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I present lattice Monte Carlo calculations for a universal four-component Fermi gas confined to 
a finite box and to a harmonic trap in one spatial dimension. I obtain the values £id = 0.370(4) 
and £id = 0.372(1), respectively, for the Bertsch parameter, a nonperturbative universal constant 
defined as the (square of the) energy of the untrapped (trapped) system measured in units of the free 
gas energy. The Bertsch parameter obtained for the one-dimensional system is consistent to within 
~ 1% uncertainties with the most recent numerical and experimental estimates of the analogous 
Bertsch parameter for a three-dimensional spin-1/2 Fermi gas at unitarity. The finding suggests 
the intriguing possibility that there exists a universality between two conformal theories in different 
dimensions. To lend support to this study, I also compute ground state energies for four and five 
fermions confined to a harmonic trap and demonstrate the restoration of a virial theorem in the 
continuum limit. The continuum few-body energies obtained are consistent with exact analytical 
calculations to within ~ 1.0% and ~ 0.3% statistical uncertainties, respectively. 
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Universal Fermi gases have been the subject of in- 
tense study in recent years. Perhaps the most inter- 
esting example of such a system in three spatial dimen- 
sions is the unitary Fermi gas: a dilute mixture of spin- 
1/2 fermions with an attractive short-range two-body 
interaction tuned to infinite scattering length. Unitary 
fermions have been realized in ultra-cold atom experi- 
ments by exploiting properties of Feshbach resonances, 
where the scattering length may be tuned arbitrarily with 
the application of an external magnetic field [TH5] (also 
see [B]). In the unitary limit, the details of the interpar- 
ticle potential become irrelevant, and consequently phys- 
ical quantities may be characterized by a single dimen- 
sionful parameter, the density of the system. 

From a theoretical standpoint, the unitary Fermi gas is 
described by a nonrelativistic conformal field theory [7J, 
and despite its simplicity, exhibits a rich and fascinating 
variety of physical phenomena. By dimensional analysis, 
the energy of the untrapped gas at density p is given by 

E(p) = £ sd E (p) , (1) 

where £3^ is a nonperturbative universal constant known 
as the Bertsch parameter [5] and E (p) is the energy 
of the noninteracting Fermi gas evaluated at the same 
density. Density-functional theory predicts that when 
confined to a harmonic trap, the energy of the system is 
given by 

E osc {Q) = V&dES sc (Q) , (2) 

where Eq sc (Q) is the energy of the noninteracting Fermi 
gas confined to a trap and Q is the total number of 
fermions, taken to be asymptotically large [5]. This re- 
sult, including subleading corrections in l/Q, has been 
been reaffirmed using a general coordinate-invariant ef- 
fective field theory description of the system [10] . 

In recent years, considerable effort has been devoted to 
determining the Bertsch parameter to high precision |llj , 



motivated in part by the fact that a theoretical determi- 
nation of may be directly compared with experimen- 
tal data (or vice verse). The current best determination 
of the Bertsch parameter from exact Monte Carlo simula- 
tions is = 0.372(5) [T2], and is in good agreement with 
the current best experimental value of £ 3( j = 0.376(4) 
|13) . Conformal and scale invariance implies an operator- 
state correspondence which relates the conformal dimen- 
sions of primary operators in free space to the spectrum 
of harmonically trapped fermions [7J . Knowledge of the 
Bertsch parameter might therefore teach us something 
about the conformal dimensions of many-body operators 
in the untrapped theory. 

Recently, Nishida and Tan [14] observed that a di- 
lute four-component Fermi gas in one spatial dimension 
is also universal when an attractive four-body interac- 
tion is tuned to produce a four-body zero-energy bound 
state. The system is conformal, and shares many prop- 
erties in common with the unitary Fermi gas in three 
dimensions [T5]. Particularly, this system is character- 
ized by a single dimensionful parameter, the density, and 
as such the untrapped system obeys the same energy re- 
lation as Eq. [l] although with a constant of proportion- 
ality £14 which need not be the same as £3^. A straight- 
forward calculation using Thomas-Fermi theory predicts 
that the trapped one-dimensional system will obey Eq. [2] 
with proportionality constant In addition to an 

operator-state correspondence, features shared by both 
the one- and three-dimensional Fermi gases include Uni- 
versal (Tan) relations, and a virial theorem for fermions 
confined to a harmonic trap [TBH213] . 

In this work, I present results from lattice Monte Carlo 
simulations of the one-dimensional Fermi gas confined to 
a finite box and a harmonic trap. I verify that trapped 
few-body systems comprising four and five fermions have 
ground state energies consistent with exact theoretical 
predictions, and also that these states obey a Viral the- 
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orem in the continuum limit. I then present numerical 
results for many-body systems, determining the value of 
Bertsch parameter £id in the infinite volume and thermo- 
dynamic limits. Although there is no known theoretical 
argument to suggest that £id and £3^ are the same, I 
show evidence that these two quantities are in fact equal 
to within 1% uncertainties. Finally, I discuss the impli- 
cations of this finding. 

The starting point for this study is a continuum ef- 
fective field theory for nonrelativistic fermions in two- 
dimensional Euclidean space-time, described by the La- 
grangian [15] : 



W-V) 4 , (3) 



where ip a is a four-component Grassmann-valued spinor 
with components labeled by cr, m is the fermion mass, fi a 
is a chemical potential coupled to the fermion number for 
each component, and g is an attractive coupling associ- 
ated with a four-body contact interaction. The contin- 
uum theory may be discretized on an N T x (2N S + 1) 
lattice with sites labeled by the integer coordinate pair 
n = (n T ,n s ), where n T G [0,N T ) and n s E [—N S ,N S ] 
following the approach of [5T]. The lattice discretization 
is carried out by identifying the various terms in Eq. [3] 
with their discrete counterparts: 



d T tp — flip 

-v 2 v> 



K 
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where b T and b s are the temporal and spatial lattice spac- 
ings, and e T and e s are corresponding unit basis vectors. 
I impose antiperiodic boundary conditions in the time di- 
rection with temporal extent (i.e., inverse temperature) 
(3 = b T N T , and open boundary condition is the space di- 
rection with spatial extent L — b s (2N s + 1). Throughout 
I work in units where H = 1. 



An arbitrary external potential may be introduced in 
the lattice construction by promoting the chemical po- 
tential to one that is space-dependent: fj, — > fj, + v n . In 
this study, I consider potentials of the form v n = (i.e., 
untrapped) and v n = k(& s 71 s ) 2 /2 (i.e., trapped), where k 
is a spring constant controlling the strength of the har- 
monic potential. In the presence of a harmonic potential 
and at infinite volume, the only relevant length-scale de- 
scribing the physical size of the system is the character- 
istic trap size Lq = l/fmc) 1 / 4 . To unify the discussion 
for untrapped and trapped fermions, I define the charac- 
teristic length scale Lq = AL/tt in the former case. 

Unitarity is achieved by tuning the coupling g to some 
critical value g c , corresponding to a zero-energy four- 
body bound state at infinite volume (with k = 0). The 
coupling may be determined by exact diagonalization of 
the lattice Hamiltonian, and yields the integral equation: 



n 



. (5) 



where A p = 2sin 2 (p/2), g c — b T g c /bg, and m = mb 2 s /b T . 
In this limit, the only dimensionful parameters in the 
theory are the lattice spacings b T and 6 S , system size Lq 
and the fermion mass m. All energy scales must therefore 
be proportional to ui = 1/(toLq) and all discretization 
errors must be a function of the dimensionless ratios e s = 
b s /L and e T = b T ui ~ e 2 . for fixed rh. 

Discretization errors in observables may be understood 
from the viewpoint of a Symanzik action [221 123j : a con- 
tinuum action in which lattice artifacts are quantified by 
couplings associated with the untuned irrelevant opera- 
tors that are allowed by the symmetries of the underlying 
lattice theory. Volume scaling may be inferred from the 
scaling dimensions of those untuned operators; following 
[24] . one expects operators of scaling dimension A to in- 
duce volume dependence in dimensionless observables O 
that scale as Lq~ a in one spatial dimension. The leading 
scaling is expected to behave as |25j : 

0(e s ) = O cont + O x t s + Oi. 666 e]- 666 + • . . , (6) 
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where the lowest order correction to the continuum re- 
sult, O con t, arises from an untuned dimension two "ef- 
fective range" term in the Symanzik action, and the sub- 
leading contribution from an odd-parity five-fermion op- 
erator with a scaling dimension A/2 = 2.333 |15j . 

For this study, I consider zero-temperature physics in 
the canonical ensemble, with fermion number q a for each 
species and total fermion number Q — ^2 a qa- From 
the fugacity expansion, one may show that the canoni- 
cal ensemble partition function Z{q) = e~^ F ^ is related 
to the grand-canonical partition function via a Fourier 
transform with respect to an imaginary chemical poten- 
tial. At the critical coupling and at sufficiently low tem- 
perature, one may define the total energy of the system 
as the total derivative: 



E( q) = hm 

,3-s.oo GtlogW 



(7) 



fixed physics 



where all physical length scales are held fixed It 
follows that the total energy for the ground-state of the 
system may be expressed as the sum of three contribu- 
tions, E = T + V + I, where 



T(q) 



v dF(q) 

hm ^rz , — 

,3-s.oo ologl/m 



/3^co O log K 



and 



dlogg c 
<91ogm 



x hm — 

/3^oo d log g c 



(8) 



(9) 



(10) 



are the expectation values of the kinetic, potential and 
interaction energy operators, respectively. In the contin- 
uum limit, it was shown that the trapped universal Fermi 
gas obeys the virial theorem: E — 2V [T5]; one therefore 
has two independent measures of the energy of the sys- 
tem, namely E — 2V and E = 2(T + I), the average of 
which yields Eq. [JJ 

Numerical simulations of the few- and many-body 
Fermi gas where performed using a fermion world-line 
representation of the canonical partition function de- 
rived in |27j . In this representation, Z(q) is expressed 
as a path-integral over all possible non-intersecting, self- 
avoiding, closed fermion loops, with the fermion number 
for each species fixed by a constraint on the winding num- 
bers in the time-like direction. Explicit expressions for 
the action and observables are provided in [27] . Metropo- 
lis Monte Carlo updating of the configuration space was 
performed using a local loop-updating scheme which pre- 
serves the constraints on the fermion paths. Simulations 
were performed using a fixed mass of m = 1.3 corre- 
sponding to a critical coupling g c sa 3.7237, and multiple 
e s and Q values to enable continuum and thermodynamic 
limit extrapolations. For each ensemble, approximately 



TABLE I. Exact and extrapolated continuum few-body ener- 
gies for fermions confined to a harmonic trap. 



Q 


E/lo (exact) E/ui 


2V/oj 


2(T + /)/w 


4 


1 1.008(6) 


1.005(10) 


1.007(11) 


5 


2.333 2.339(7) 


2.331(12) 


2.346(13) 



150-350 uncorrelated configurations were generated after 
thermalization. 

From properties of the Schroedinger algebra one may 
show that the spectrum of few- and many-body systems 
confined to a harmonic trap consist of a tower of levels 
separated by 2u> [7 . Trapped simulations were therefore 
performed at a temperature (3uj > 10 to ensure adequate 
suppression of excited state contamination. In the case of 
untrapped many-body studies, the energy splittings are 
expected to be of order the Fermi energy, E F = Q 2 ui/2, 
and therefore for each to the temperature was chosen such 
that /3E F > 10. 

For four and five fermions confined to a trap, esti- 
mates of the energies in units of w are plotted in Fig. [I] 
as a function of the discretization error e s . Continuum 
extrapolated energies obtained from a three-parameter 
least-squares fit to Eq. [6] are provided in Table |T] and 
are consistent with exactly determined values of unity 
and 2.333 to within 1.0% and 0.3% uncertainties, respec- 
tively. All fits yielded a \ 2 P er degree of freedom (d.o.f.) 
of < 1. The exact continuum energies were determined 
via the operator-state correspondence using scaling di- 
mensions for few-body operators calculated in [T5, . Con- 
vergence of the energies defined by 2V and 2(T+ 1) in the 
continuum limit demonstrates a restoration of the virial 
theorem. 

The one-dimensional Bertsch parameter is defined as 
the ordered limit £i<j = limQ^oo lim es _>o £q( £ s) 7 where 
£q(cs) is the (square of the) energy of the untrapped 
(trapped) Q-body system at finite e s measured in units 
of the free-gas energy, E = Q 3 lo/2A {E° sc = Q 2 w/8). 
Untrapped simulations were performed at twelve equally 
spaced values of k F bs = Q( s € [0.15,0.7], where k F is 
the Fermi momentum, and for fermion numbers Q = 
32,48,56,64,72,80,88. Trapped simulations were per- 
formed for all integer values of l/e s G [7,20] and total 
fermion number Q = 28,32,36,40,44,48,52,56 subject 
to the constraint Q 1 / 2 e s < 1.0. Continuum extrapola- 
tions of the energy defined by Eq. [7j were performed at 
fixed Q over those respective intervals using Eq. [6] with 
three fit parameters. A x 1 /d.o.f . < 1.5 was achieved over 
each fit interval and variation of the maximum e s value 
used in the fit yielded statistically consistent extrapo- 
lated values for the energies, indicating that the fits were 
robust. Thermodynamic limit extrapolations of £q(0) 
were subsequently performed. For the untrapped system, 
an extrapolation in 1/Q was performed using a constant 
plus linear ansatz fit function, yielding t;id = 0.370(4) 
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FIG. 2. Thermodynamic limit extrapolation of the Bertsch 
parameter for untrapped and trapped fermions. 



with a x 2 /d.o.f. ~ 1.6. For the trapped system, an ex- 
trapolation was performed using a constant fit function, 
yielding f ld = 0.372(1) with a \ 2 /d.o.f. « 0.3. Contin- 
uum limit extrapolation results for £g(e s ) are shown as 
a function of 1/Q in Fig.[2j along with fit results and er- 
ror bands for the thermodynamic limit extrapolation of 
£q(0). The values obtained for £i<j are statistically equal, 
and moreover, are consistent with the best known Monte 
Carlo, and experimental values for the three-dimensional 
Bertsch parameter £3^. 

Provided the interparticle potential is attractive and 
the system is stable, one may argue that the one- and 
three-dimensional Bertsch parameters must lie between 
zero and unity. The likelihood that the agreement dis- 
covered between £1^ and £3^ is purely a coincidence is 
roughly 1/100. From the point of view of the operator- 
state correspondence, this equality seems to imply a sur- 
prising connection between the conformal dimensions of 
many-body operators in the one- and three- dimensional 
universal Fermi gas. Whether this curious relationship is 
due to symmetries or is dynamical in origin remains an 
open question worth pursuing. If this equality is indeed 
exact, then numerical simulations of the one-dimensional 
Fermi gas could offer a computationally inexpensive al- 
ternative for computing £3^ to high precision. Further 
exploration of the one-dimensional Fermi gas is under- 
way, including few- and many-body studies of the inte- 
grated contact density, and may provide deeper insights 
into this intriguing relationship. 
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